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Impact of membrane bistability on dynamical response of neuronal populations 


Wei Wei^’^, Fred Wolf^ and Xiao-Jing Wang^’^d 
^ Center for Neural Science, New York University, New York, NY 10003 
^ Department of Neurobiology and Kavli Institute for Neuroscience, 

Yale University School of Medicine, New Haven CT 06520 
® Max Planck Institute for Dynamics and Self-Organization and Bernstein 
Center for Computational Neuroscience, D-31011, Cottingen, Germany 
NYU-ECNU Institute of Brain and Cognitive Science, NYU Shanghai, Shanghai, China 

Neurons in many brain areas can develop pronounced depolarized state of membrane potential 
(up state) in addition to the normal hyperpolarized down state near the resting potential. The 
influence of the up state on signal encoding, however, is not well investigated. Here we construct a 
one-dimensional bistable neuron model and calculate the linear response to noisy oscillatory inputs 
analytically. We find that with the appearance of an up state, the transmission function is enhanced 
by the emergence of a local maximum at some optimal frequency and the phase lag relative to the 
input signal is reduced. We characterize the dependence of the enhancement of frequency response 
on intrinsic dynamics and on the occupancy of the up state. 
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Introduction- Elevated state of neuronal membrane po¬ 
tential (MP), the so called up stafe has been observed 
extensively in different brain areas [l| . In this regime the 
MPs of neurons are characterized by a bimodal distribu¬ 
tion resulting from two stable fixed points in membrane 
dynamics. This bistability of neuronal dynamics leads to 
synchronous transitions between the down and up states 
of neurons in a network and the development of global up 
and down states ii- The exact role of the bistability 
of membrane potential in signal encoding and processing 
is still not well understood. 

Individual neurons in a network receive noisy synaptic 
inputs and fire spikes irregularly [^. Besides the sta¬ 
tionary firing rate, one important characteristic of neu¬ 
ronal dynamics is the response to time varying signals 
superimposed on background noise Q. This dynamical 
response of cortical neurons has recently been measured 
experimentally up to 1 kHz of signal frequency, which 
revealed very hight cutoff frequencies 0. Theoretically, 
the linear response has been obtained analytically for the 
leaky integrate-and-fire (LIF) neuron Q, in which mem¬ 
brane dynamics has only one stable fixed point, and for 
the r — T model (a piecewise linear version of experien¬ 
tial integrate-and-fire model) Q, in which an additional 
unstable fixed point for action potential initiation was 
included. The effect of the unstable fixed point in mem¬ 
brane dynamics on dynamical response was also investi¬ 
gated numerically in other one-dimensional models and 
conductance-based models Q. Theoretical characteriza¬ 
tion of dynamical response of neurons that exhibit up and 
down states, however, is still missing. Intuitively, when 
the membrane potential of a neuron has a higher proba¬ 
bility to be around some depolarized voltage, it is more 
likely that a small oscillatory signal can contribute to the 
firing of spikes, leading to enhancement of frequency re¬ 
sponse. Here we propose an analytically solvable bistable 


neuron model to investigate the impact of the up state 
on the dynamical response of neurons. 

Model description- In this work we construct a one¬ 
dimensional bistable neuron model, which has piecewise 
linear subthreshold dynamics and is analytically solvable 
for the linear response. The dynamics is described by the 
following Langevin equation, 

TV = f{v) -\- fj.-\-ar]{t) , (1) 

where 

{ — V , —OO < V < Vo 

ri{v - vti) , vq <v <vi ( 2 ) 

r{v - vto) , vi <v <Vb 

Here v is the MP relative to the resting potential, p- is 
the mean external input, ri{t) is a Gaussian white noise 
which satisfies {r]{t)) = 0 and {ri{t)ri{t )) = T6{t — t ), 
and a is the strength of the noise, r is the membrane 
time constant near the resting potential, which we will 
take as the unit of time in the theoretical result. Note 
that the membrane dynamics here might result from an 
interaction between patterned synaptic inputs and intrin¬ 
sic membrane dynamics. Fig. lA shows an illustration 
of the model dynamics when there is no external input. 
The slopes of the middle and right pieces are denoted 
as ri and r, respectively. Note that r characterizes the 
membrane dynamics around the higher stable fixed point 
where the time constant is given by T/|r|. MPs at the 
crossing points of the left piece with the middle piece, 
and the middle piece with the right piece are denoted as 
Vo and vi, respectively. When there is no external input, 
Vti and Vto are the unstable fixed point and the higher 
stable fixed point, respectively. We will fix vq and Vto, 
whereas Vti and vi are given by Vti = (1 -f l/ri)uo and 
vi = (rivti — rvto)/(ri —v). The deterministic dynamics 
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FIG. 1: Illustration of the model. A, illustration of the piece- 
wise linear model; B, MP trajectories for ri = 1, 5 and 10 
from bottom to top. Parameters used are: r = —1, vo = 
0.5, Vr = Vti, Vto = 2, Vb = -0.2, r = 10 ms, Tr = 
0 ms, = 0, a = 0.5. 


(cr = 0 in Eq. ([T])) possess one lower stable fixed point, 
one unstable fixed point and one higher stable fixed point, 
located at v = fi, Vti — /i/ri and Vto — respectively. 
When the MP reaches an absorbing boundary Vb, it is 
reset to a resetting potential Vr for a refractory period 
Tr- The larger one between jilol and l-Dbl determines the 
rheobase current of the model neuron, where Vg = —Vq 
and Vb = r{vb — Vto). Since the model neuron will never 
fire spikes when |r| is very large if Vb is fixed, we will fix 
Vb and choose Vb determined by Vb = Vto + Vb/r. The 
relative values of vg and Vb might indicate different dy¬ 
namical regime of the model neuron. For example, a 
larger I'Obl implies that the neuron can make a transition 
from the up-down regime to a tonically depolarized state 
when fx > |wo|, which is reminiscent to that observed 
for cortical neurons from sleep to wakefulness [l^. We 
are most interested in the noise-driven regime, i.e., the 
mean external current is smaller than the rheobase cur¬ 
rent, since real neurons work in a regime in which excita¬ 
tory and inhibitory synaptic inputs to individual neurons 
balance each other 11| . In this dynamic regime, our neu¬ 
ron model describes barrier penetration in a double well 
with reinjection of probability current after reaching an 
absorbing boundary. Here we focus on the case when Vr 
is located within the middle piece, therefore the model 
neuron could fire several spikes before transiting to the 
down state, as normally observed for neurons exhibiting 
up state. Fig. IB shows the MP trajectories for three 
different ri. We see that with a larger ri, the MP spends 
more times in the up state. 

The Fokker-Planck equation (FPE) framework- The 
FPE corresponding to Eq. m has the following form 


dtP{v, t) + dy{f{v) F n - Ddy)P{v, t) = 0 , (3) 

where D = is the diffusion constant. We will use 
both D and cr in the following. Defining the probability 
current as J{v,t) = {f{v) + n — Ddy)P{v,t), the FPE 
then becomes the equation for probability conservation, 
dtP{v,t) + dvJ{v,t) = 0. 

The boundary conditions are specified in the follow¬ 
ing (subscripts 1, 2, 3 indicate the left, right, and mid¬ 


dle MP regions in Fig. lA). At the absorbing boundary 
Vb, P 2 {vb,t) = 0. At the resetting point Vr, P 3 {Vy,t) — 
P?,{vf,t) = t) aaddyP3{v+ ,t)-dyP3{yf ,t) = dyP2{vb,t- 
Ty), from the the resetting condition and continuity of 
the probability density and probability current. At no 
and vi , the probability density and its derivative are con¬ 
tinuous: Pi{vQ,t) = P3{vo,t), dyPi{vo,t) = dyP3{vo,t), 
P 3 ivi,t) = P 2 (vi,t), and dvP 3 {vi,t) = dvP 2 {vi,t). Fi¬ 
nally the normalization condition of the probability den¬ 
sity requires lim.u_)._oo Ti(n, f) = 0. With these bound¬ 
ary conditions the asymptotic solution of the FPE is 
uniquely determined (the possible transient is not of in¬ 
terest here). The instantaneous firing rate is given by 
the probability current through the absorbing boundary, 
v{t) = J{vb,t) = -DdyP2{vb,t). 

Development of the up state- When the mean input 
to a model neuron is constant, the stationary probabil¬ 
ity density, denoted as Po{v), can be obtained by setting 
J(v,t) = Vg 0 . Here vg is the stationary firing rate, 
which is determined by the normalization condition of the 
stationary density, fl'’^Po{v) dv = 1. The existence of 
the up state requires the appearance of a local maximum 
of the probability density at a depolarized MP value. 
Two peaks appear in Po('*^) if there exists an up state 
in the MP trajectories, located around the lower stable 
fixed point and the higher stable fixed point, which will 
be denoted as and Pq^, respectively. The MPs 

corresponding to the two peaks, denoted as Vdown and 
Vup, are the mean values of MPs at the down state and 
up state, respectively. From the expression of Pg{v) 0 : 
it is easily to see that the down state locates at the lower 
stable fixed point, Vdown = The development of a local 
maximum around the higher stable fixed point requires 
Pmi'^) = 0 since the probability density decreases mono- 
tonically when 0 < f < ui, or equivalently, the following 
equation 


xe 




( 4 ) 


has a solution within the range vi < v < Vb, where x = 
r{v^)+u ^ ^ith Vup independent of ri. 

Since Xb < x, we have x > 0 and vi < Vup < vtg — pjv 
from Eq. Q. Therefore the mean value of the MP at 
the up state locates lower than the higher stable fixed 
point in the deterministic dynamics due to the influence 
of noise and the absorbing boundary. The probability 
density at v = Vup is given by 


pup 


r{vup - Vto) + h- ' 


( 5 ) 


The dependence of Po{v) on ri and r is shown in Fig. 
2A-B. With the increase of ri, the transition from the up 
state to the down state becomes more difficult, therefore 
the MP resides on the up state for a longer time (Fig. 

indicating a larger ratio between the maximal prob¬ 
ability density at the up state and down state (Fig. 2A). 
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Note that changing the slope ri can determine whether 
the up state exists or not by adjusting ui, but does not 
influence the position of v^p if it exists. With the increase 
of |r| (the absolute value of r), v^p is shifted slightly to¬ 
wards the higher deterministic fixed point and the ratio 
pup/pdown decreases (Fig. 2B). 

Linear response- Now consider a weak sinusoidal signal 
encoded in the mean input, pL{t) = pLe cos{ujt), where £ 
is small. At the linear order in e, the instantaneous firing 
rate is given by v{t) = vq e:|:/ic(a;)| cos(wt — (pc^uj)), 
where |z/ic(a;)| is the transmission function and (pd^) is 
the phase lag. We find that a complex response function 
iyic{u!) can be obtained analytically by solving the FPE 
at the linear order using the Green’s function method. 
The transmission function is the absolute value of vidoj), 
while the phase lag (pd^) is given by the phase angle, 
(/)c(w) = arg(z^ic(w)). The expression of i^ic(w) reads 


i/ic(w) = 


B 


iuj{\ -I- 1/ri) 

(1 — za;)(l -b iuilri) 


{diPai - V^$iFoi) 


iui(l/ri — 1/r) 

(1 -b iuj/r){l -b iiojri) 
d Dr I (r - ri) 


(■ 01^5 - diYdPo2{vi) e^^ 


(ri + iu}){r -b iuj) 
Vo/d Dr I 


- i>iP5i)Po2di) e 


Ai 


, (V’ln'lr - d'lY,,r)e^° - X 

1 + tuj/ri I luj/r 


ddiY, - diYdY^ - idiYo - diYdY2)e^ 


( 6 ) 


where tpid), ‘hi('y), etc. are parabolic cylinder functions 
0 , and Fi(r;), Yo(v), B, etc. are combinations of them 
to simplify the expressions, as defined in 1^. Note the 
functions adopt their values at v = vq, unless denoted 
otherwise. Taking w —>■ oo in Eq. dSD, we find that 
the high frequency limit is the same as the LIF model 
and the r — t model, i.e., vic —>■ high 

frequency limit is characteristic of the linear dynamics 
and absorbing boundary Q. 

Frequency-selective enhancement of linear response by 
the up state- The dependence of the linear response on 
ri and r are shown in Fig. 2C-D. We see that a local 
maximum, denoted as appears at frequency 

(Fig. 2C-D, upper panels) accompanying the develop¬ 
ment of the up state (Fig. 2A-B). The local maximal 
value of transmission function at the resonance frequency 
increases with ri and decreases with |r|, following the 
same trend as the ratio between the probability density 
at Vup and Vdown (Fig. 2A-B). The phase lag of the fir¬ 
ing rate response relative to the input signal is reduced 
when there is a more pronounced up state (Fig. 2C-D, 
lower panels). Therefore the up state can enhance the dy¬ 
namical response by developing local maximum at some 
specific resonance frequency, and reduce the phase lag of 
the response. 

We characterize the relationship between the up state 
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FIG. 2: Dependence of probability density and linear response 
on n and r. A and B, probability density for different n and 
r. C and D, dependence of the transmission function (upper 
panels) and phase lag (lower panels) of the linear response on 
ri and r. The linear response is normalized with the value 
at / = 1 Hz. Signal frequency / is related to the angular 
frequency w by w = Stt/. Solid lines are from theoretical 
results and asterisks are from simulations. Parameters used: 
r = —1 in A and C, ri = 10 in B and D. Other parameters 
are the same as in Fig. [T] 


occupancy and resonance in the linear response quanti¬ 
tatively in Fig. 3. While the ratio between probability 
densities at v^p and Vdown increases with ri (Fig. 3A), 
the maximal value of the transmission function at the res¬ 
onance frequency also increases with ri (Fig. 3C). This 
leads to an increase of with p^P/pdown ('p-g_ 
inset). Similarly, v^d^ decreases with |r| (Fig. 3D), fol¬ 
lowing with the same trend as p^P /pdown g^cept for an 
initial small |r| regime where there is no local maximum 
for |izic(w)| (Fig. 3B). This also leads to an increase of 
udd^ with p^P/pdo^^ (Fig. 3D, inset). The signal fre¬ 
quency at which the transmission function is maximally 
enhanced, keeps constant when ri increases, i.e., 

being independent of the time scale characterizing the un¬ 
stable fixed point. On the contrary, goes to higher 

frequencies when |r| increases, therefore is determined by 
the membrane time constant at the up state, T/lrj. 

We tested the predication by using a neuron model 
with biophysically realistic membrane currents, which in¬ 
cludes a noninactivating potassium current controlling 
the up state and an inward rectifying potassium current 
stabilizing the down state Q. We find that the normal- 
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FIG. 3: Dependence of up state occupancy and resonance 
in the linear response on ri and r. A and B, the up state 
Occupancy p^P ^as a function of n (A) and |r| (B). 
C and D, the local maximal value of the normalized 

transmission function as a function of ri (C) and |r| (D). E 
and F, the resonance frequency as a function of ri (E) 

and |r| (F). Insets: as a function of P^p/P§°'^^ due 

to the change of ri (C) and and |r| (D). Parameters used: 
r = — 1 in A and ri = 10 in B. Other parameters are the 
same as in Fig. [T] 


ized frequency response is strongly enhanced with the 
development of the up state [14 1. 


Discussion- Dynamical response to noisy oscillatory 
inputs is a fundamental characterization of neurons’ 
capability in signal encoding and transmission. Our 
model predicts a frequency-selective enhancement of sig¬ 
nal transmission with the development of the up state. 
This prediction can be experimentally tested for neu¬ 
rons exhibiting up and down states by performing ex¬ 
periments similar to that in Q. Selective enhancement 
of specific high frequency component of signals might be 
functionally beneficial for neural communication, since 
signals with frequencies within the gamma band (30- 
100 Hz) or “high gamma” (> 100 Hz) have been sug - 
gested to synchronize inter-regional brain activity [l5|. 
Bistable piecewise linear membrane dynamics was intro¬ 
duced previously to approximate the nullcline of MP in 
the FitzHugh-Nagumo model [l^. In the limit of r —>■ 0, 
the linear response of that system was obtained [3. 
Here we obtain the linear response analytically for a one- 
dimensional bistable system with general r. Bistability 


in individual neurons provides an alternative explana¬ 
tion for persistent activity in some brain areas (l8j |. and 
is related to perfect temporal information accumulation 
0 and self-organized criticality [i^. Our work sheds 
insights into the possible role of the up state on signal 
encoding and transmission through population response. 
Further work on building a network of interconnected 
such bistable units is needed to examine the relationship 
between the up state at the individual neuron level and 
at the circuit level [2lj |. 

We thank J. Mejias, J. Murray and F. Song for care¬ 
fully reading the manuscript. 
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